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ABSTRACT 


Adaptive control system techniques have been utilized to inves 
tigate identification of the dynamic equations of operating control 
systems. The identification has been achieved by using model modu 
lation techniques to drive dynamic models into correspondence with 
operating control systems. The system identification then proceeded 


from examination of the model and the adaptive loop. 


The model modulation techniques applied to adaptive control 
systems are breifly discussed These techniques applied to simulation 
studies on an analog computer and a digital computer, and applied 
directly to a position servo-mechanism, are then discussed The 
discussion presents data showing the accuracy of identification of 
the dominant characteristics of the dynamic equations when the order 
of the equation is either known or approximated. It was found that 
these dominant roots could be found within an average accuracy of 5% 


for both complex phase angle and undamped natural frequency. 


The authors wish to express their appreciation to Professor 
Richard C. Dorf of the U. S. Naval Postgraduate School, for the 
suggestion of the model modulation method and for his assistance 


during this investigation. 
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CONTROL SYSTEM IDENTIFICATION THROUGH MODEL MODULATION METHODS 


Chapter I 


INTRODUCTION 


General 

The last three decades have seen control system technology grow 
from a series of rules of application to a full regime encompassing 
most of the disciplines of engineering. The study of control systems 
has progressed deeply into analysis and synthesis of linear and non: 
linear control systems. There are, additionally, offshoots of basic 
Systems with most exhaustive methodology: examples are the sampled 


data, hybrid analog-digital and multi-loop self adaptive systems. 


The problem of control system identification led to some of 
the original investigation in the control field. Despite the increase 
of knowledge in the field of control engineering, system identification 


has remained amajor problem. 


The intent of this project was to investigate a method of control 
system identification utilizing techniques from the growing technology 
of the self-adaptive control system. One restriction was laid upon 
the project at its inception. This was that the method of system 
identification remain independent of the system to be identified 
Ideally then, an operating system should be identified solely through 
examination of the input and the response. To this end a control 
system identifier has been modeled and tested on an analog computer, 


through digital computer analysis, and upon a position control System. 





Control System Identification 


The characteristics and performances of the fundamental components 
of feedback control systems are generally amenable to systematic analysis 
The control system, when composed of these fundamental components, may 
be easily synthesized and analyzed. Components, however, when removed 
from the theoretical state do not exhibit the invariant characteristics 
attributed to them. Furthermore, control systems are often synthesized 
on the basis of transfer function manipulation. This synthesis occurs 
under the assumption of negligible loading of one component system by 
another. This assumption is too often fallacious. In either or both 
of these situations the control system engineer is confronted with a 
control system which does not conform to previous analysis. And if 
this analyzed control system is a portion of some larger process, it 
may be necessary to reevaluate the control characteristics and perform- 


ance, in other words to identify the system. 


In some extremely variant systems the processes of system identi- 
fication are utilized to serve as controllers in a self-adaptive control 
System. The most notable of these cases are the self-adaptive control 
Systems currently in use and being further developed for aircraft 
installation. The NASA X-15 is an example. In this installation the 
response characteristics are examined and the results are used to vary 
flight control loop gain to maintain near invariant control response. 
The equipment used for this process is quite complex, but is justified 


by the operating environment of the vehicle. 





For less complex control systems or in less complex environments, 
the problem of control can be solved by more conventional means, without 
recourse to dynamic self-adaptive system technology. These less complex 
systems require accurate identification of the process dynamics. Once 
the transfer function of a process has been evaluated this evaluation 
can serve as a basis for design of the remainder of a system or for 


readjustment of the process through compensation techniques. 


Identification Methods 

Truxal writing in reference (1) has placed the methods of control 
system identification into four broad categories: 

a. Parameter identification. 

b. Frequency response methods. 

c. Impulse response methods. 


d. Direct evaluation of system differential equations 
Of motion. 


Some of the difficulties in one-by-one determination of model 
parameters have been mentioned. Additionally, this method is possible 
only when the detailed operation of a process is reducible to funda- 
mental differential equations in which the primary parameters may be 
included and the negligible ones discarded. This ideal situation 


rarely exists for complex control of complex processes. 


As frequency response techniques are well-known, a description 
of them has been omitted. These procedures are extremely important 
in general control design. This is the most common method of iden- 
tifying the transmittive dynamics of a control process. Many working 
control systems however, cannot be subjected to variable frequency 
oscillations for reasons of system complexity, operating economy, or 


— 





because of the extensive time and effort required for determination 
of the gain and phase characteristics. Two additional problems occur. 
First, in some circumstances, frequency response testing should not 
be superimposed on regular operation, because testing would interfere 
with designed operation. Secondly, a working system is subjected not 
only to its driving input, but also to noise occuring during the op- 
erating process. The extraneous disturbances of normal operation 


would make the interpretation of frequency test results difficult. 


The impulse response can serve to completely identify any linear 
time-invariant two port system. The various impulse response techniques 
are utilized in many of the more sophisticated adaptive control systems. 
Truxal has cited three primary techniques for measuring impulse response. 
The first and most direct method is to excite a system with an impulse 


and measure the response. Second, the system can be excited by a 


known signal. For an initially inert system, the impulse response 
З (є) can be determined from the input, ms and output, 8 (c)? through 
the convolution integral 
Т 
(t) Jt) * (t- T) 2! 


° 
Finally the impulse response can be measured by application of a small 


random signal to the process. If white noise is used as the input, 
the cross correlation function between the output and the noise is the 


impulse response. 


The attempt to utilize the impulse response in identification is 


complicated by three problems. Noise has a much more delerterious 





effect on the output than in the case of frequency response tests.  Non- 
linearities are much more difficult to detect than in frequency response 
tests. Finally the mechanization necessary to determine the impulse 
response is almost prohibitive for use in system identification. When 
the impulse response method is used in adaptive system control the 
engineer is generally interested not in the full identification, but 

in identifying only some characteristic of the response. Quantitative 
measurements are made to serve as sufficient identification to allow 
adaptive correction. As an example, the number of oscillations during 

a specific period of a system response would identify the impulse 
response. These approximations are not sufficient for full identifi- 


cation, however. 


The evaluation of the differential equation for a linear process 
offers an additional method of identification. The differential equation 


of a process might be expressed as an equation of input 2 and output 


S CB in the following form: 
N ሎሠ/ ክላ 
de. + a, e. 4 seo й, ФС + A, = kaar Б 1.22 


The a, term may be evaluated from the nature of the zero frequency or 
steady state gain. If the order of the differential equation is known, 
and the nature of the right hand side is known, the other 84-1 po 
efficients may be evaluated by generalizations such as minimizations 
of square of the error e, where e is defined: 


de, LO 4 2,46 * q C B ኪና in, 


(9 
ዘ 





The ane] coefficients are evaluated through use of this error criterion 
by generation of successive differentiations of the response. Unfor- 
tunately, searching over an (N-1) space for a minimum is a problem of 


potential difficulty. General problems also arise in settling at 


local minimums and in error system stability. 


Many of the problems of systematic determination of the differ- 
ential equation coefficients are alleviated by the use of system-model 
reference. This is a natural variant to the adaptive control system 
technology. The exploitation of this method involves utilization of 
the difference between system and model response to the forcing input. 
This difference has been handled by several methods in present technology. 
In general self adaptive processes, the difference as operated upon 
by an adaptive technique is used to drive the system into operating 


correspondence with the invariant response of the model. 


This investigation was based on a variation of this method. During 
this investigation adaptive techniques have been used to drive the 
model into operating correspondence with the operating system. The 
characteristics of the model differential equation have been modified 
by the adaptive loop. These characteristics then have been used to 


accurately define the coefficients of the system differential equation. 


Chapter II 


MODEL REFERENCE ADAPTIVE SYSTEMS 


Introduction 

A major advantage of the model reference adaptive system is that 
direct measurement of the operating system differential equation 
parameters is not required. This factor results in considerable 
economy in mechanization of the adaptive components. As previously 
mentioned there are in existence various methods of self adaptive 
control. However, in many of these systems the adaptive loop mech- 
anization requires the assumption of the second order nature of the 
system. An example of this is Osder's system detailed in reference (1). 
For this system, the sign reversals in the output oscillations are 
counted and the count used to correct the second order zeta. Most of 
the impulse response self adaptive systems also require assumption of 
a basically second order system to keep the equipment within allowable 
limits. The model reference system relieves this restriction. The 
model reference adaptive system can be utilized for higher order systems 
about which the second order assumption is invalid. Additionally, 
current literature presents several examples in which model reference 


allows the control of several parameters within the operating system. 


References (1) and (2) present several methods of model reference 
adaptive control systems which have been proposed or developed. The 
variants of these systems include the Ham-Lang conditional control 


System, which in analysis is passive adaptation through model reference. 
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The other extreme is perhaps, Braun's method, which utilizes an ortho- 


normal spectrum analyzer to solve the identification problem. 


Parameter Perturbation 

Perhaps the most straightforward adaptive mechanization results 
from the use of model reference with perturbation techniques. McGrath, 
Ra jaraman and Rideout in reference (3) have proposed and examined 
various parameter perturbation model reference schemes. Figure 1, 
taken from their discussion presents two variants of the parameter 
perturbation scheme. 


Fig. 1 
PARAMETER PERTURBATION METHODS 


MODEL 





The adaptation processes in these systems are driven by the error, є, 


between the model and system. This error is used to form a performance 


— 





criterion, which for this example was integral of error squared, Б La 


This criteria contaíns a component of the frequency at which a parameter 
of the system in Fig. l(a) is being perturbed. The amplitude and 

phase of this component give the magnitude and sign of the signal 
obtained through perturbation and integration. This signal is then 

fed back negatively to reduce the error. The adaptive loop mechan- 
ization is illustrated in Fig. 2. 


Fig. 2 
PARAMETER PERTURBATION ADAPTIVE LOOP 


E A 





For analog simulation electronic multipliers have been utilized for 
Parameter correction insertion. In operating systems servo-multi- 


pliers could be utilized. 


The system of Fig. l(a) suffers from the disadvantage that the 
perturbation appears in the output. The system of Fig. l(b) avoids 
this problem through model parameter perturbation. This second method, 
however, requires that the model be nearly the same form as the system. 
Both methods, moreover, require phase correcting networks to correct 
the phase shift between parameter perturbation signal in the model and 
the oscillator signal. Through utilization of further model forms 


with parameter perturbation rhe authors demonstrate that it is possible 


=o 





to prepare a configuration which is both system and signal adaptive, 
using multiple models with perturbation of either model or system 
parameters. Rajaraman in a later article, Ref. 4, has shown the 
detailed operation of the signal and system adaptive process. In 
summary, this signal and system adaptive process requires two models 
and two adaptive loops. Тһе two models are of different nature. The 
first is an ídeal version of the process, and the second is of the 
same order as the process. The second model is optimized for the 
input by the first adaptive loop operating upon the difference 
between the output of the two models. The second adaptive loop 
causes system parameters to follow the optimized second model par- 


ameters, and at the same time corrects for disturbances. 


Model Output Modulation 


Dorf and Byers have shown the feasibility of the model reference 
self adaptive system using output modulation techniques. The analysis 
of the techniques and a discussion of experimental results are found 
in Ref. 5. They are briefly summarized here. Fig. 3 presents the 
block diagram of the system. 


Fig. 3 
MODEL MODULATION ADAPTIVE LOOP 


To PARAMETER | CoRRECTION 


Demop- 
ULATION 








In this adaptive loop, the output error is: 


N | * A sin - 


For this example error squared has been used as the criterion. Thus, 
after being shaped and demodulated a phase-sensitive d.c. signal is 
obtained. When this signal is subjected to integration and filtering, 


the open loop output of the adaptive circuit becomes: 
А С በህ. ሚገቡ 22 
m m 5 


The adaptive circuit loop was closed by this signal being inserted 
into the system as a driving function, to correct a system parameter. 
The result of this multiloop circuit is a complex, non-linear system 
which has adaptive characteristics to maintain an approximately in- 
variant system transfer function. In functional notation the system 
output might be described as a function of the following variables: 


Cry BE Ka Y Б Tes Т, М) 


(t) 
Kid = adaptive loop gain 
v. = frequency of input signal 


R = magnitude of input signal 


Te = filter network time constants 
T = system transfer function, uncorrected 
M = model transfer function 


The adaptive loop filter was designed to filter out the normal system 
operating signals. The speed of adaptation depends partially on the 
filter transfer function. The gain and lag of the filter influence 
the magnitude of the adaptive response and the input signal contami- 
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nation. One important feature of this method is that an absolutely 
constant frequency oscillator is not required. The operation of the 
adaptive loop is not dependent upon frequency stability of the modu- 
lating frequency. Furthermore, as there are no significant time 
constants between modulation and demodulation, the frequency of 
modulation і5 limited only in bounds at the lower frequency by the 
necessity of being an order above the input signal frequency. Dorf 
and Byers further investigated a two parameter correction situation, 
driving one adaptive loop with output and the second adaptive loop 
with the derivative of output. In their example the system open 


loop transfer function was: 


я 
Kw 


CXS =. a ирт 


s(s*+ 2Sun, + wi) "S 


The model open loop transfer function was similar. The first adaptive 
loop was driven by С. апа E to correct (K; + AK), the open loop gain. 
The second adaptive loop was driven by С апа C. to control the inner 
loop damping (2 G o + A 2 € ш). The authors concluded that a two 
parameter self adaptive system of this kind was feasible. A major 
restriction was found to exist if the two parameters being adapted 
were not independent, and the adaptation would be successful only if 


the range of required adaptation was restricted. 


Application of Model Modulation 


To verify the results of model modulation, and to obtain exper- 
ience with the characteristics of the method, the method was simulated 


on a Donner 3100 analog computer. Initial investigation concerned 


-12- 





the use of a second order system with one variable parameter being 
adapted to follow a second order model. The model and the system 


were ordered for an open loop transfer function: 


E 


o E .]oo . 9n 
s S(s t2 Sur S (s 10) 


The parameter drift then is represented as 


e 


G(s) = — 2. 
S(S + — +n) 


where n was a step parameter drift. The input was also of various 
wave forms of frequency less than the bandwidth. The adaptive system 


open loop is represented as: 


2 
"A. 
G(s) =  - ዚ ሔ << - ሩዙ 2» 
(сетш чу-е m.) 
where m. is the correction derived from the output driven adaptive 


loop. This system might be considered to simulate a second order 


position servo as displayed in Fig. 4. 


Fig. 4 
SECOND ORDER SERVOMECHANISM 





Constant Armature Current D.C. Motor Drive 


— 





Thus the variable drift simulated a system with varying or unknown 
friction(f), a system with constant gain but varying velocity error 
coefficient. The correction would be inserted by varying tachometer 


gain feedback magnitude. 


For operational simplicity the error criterion was altered from 
the integral of error squared, ISE, to the integral ot absolute 
error or IAE. Additionally, the modulation and demodulation was 
accomplished through use of electronic multipliers. incidental noise 
problems with these multipliers was such that the pure integration 
followed by a filter was removed from the adaptive loop, and ínstead 
a double filter network substituted. The alteration is shown symbol- 
ically in Fig. 5. 


Fie.) 
MODIFIED MODEL MODULATION LOOP 


ABSOLUTE To PARAMETER 


Cen — МИ St CORRECTION 


ነ 
&«b 








| А SIN wt 


Satisfactory results were achieved with this circuit with the filter 
break points a = .1, and b= .1. The full computer diagram is shown 

in Fig. 6. The value of the gain of the adaptive loop for the circuit 
of Fig. 6 was 3. Fig. 7 presents the output data for this same circuit, 


illustrating adaptation for various parameter drift conditions 





Additional investigation was made of the control of a third order 


system with one variable parameter by a second order model. The third 


order system investigated open loop transfer function was 


K_ 100 


G(s) * s (s + E (s + 10 + n - m.) 


The model utilized was that of the previous example. The analog 
computer setup was that of the transfer function format similar to 


Fig. 6, with an additional function K. installed. The adaptive 
s + K 
р 


loop was unchanged. This system then simulated a third order control 


System as illustrated in Fig. 8. 
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For this situation the variable parameter simulated would be p, 

with the correction being inserted through a small servo varying m. 
The adaptation was satisfactory for K, = 40, when p = 10. Adapt- 
ability became difficult for K, < 40 for large parameter shift. Fig. 9 


presents the output data for this circuit for a representative K, - 40, 


under various parameter drift conditions. 


As a final illustration of the model modulation method, cancel- 
lation techniques were used for system control. The process is best 


illustrated in Fig. 10. 
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Figure 9a 
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Fig. 10 
CONTROL SYSTEM CASCADE COMPENSATION ADAPTATION 
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This technique is most useful when the operating system pole is not 
amenable to correction because of the physical characteristics of the 
system. In the operation illustrated by Fig. 10 zero drift caused 
compensator cancellation. When drift was introduced the correction 
forced the zero of the compensator to cancel the system pole, and the 
desired system operation was restored. The analog computer operating 
schematic is presented in Fig. 11. For operational simplicity the 
drift was inserted as a step change in the system pole. The computer 
output data is presented in Fig. 12, where the successful adaptation 


of the system is well illustrated. 
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The results of these simulations of the adaptive system were 
considered to satisfactorily demonstrate the feasibility of the 
method and to serve for participant training. The analytical proof 
of the correctness of the alterations to the initial formulation, 
IAE for ISE as criterion, and filter instead of pure integration, 
have been neglected. The demonstrated accuracy of the adaptation 


was felt to be sufficient proof of the method. 
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Chapter III 


IDENTIFICATION THROUGH MODEL MODULATION 


The methods of successful adaptive control through model modu- 
lation were directly applicable to the identification ot operating 
systems. The requirements initially set were that the identification 
Should proceed utilizing only the operating system input and output 
without disturbing the system operation. To utilize the model modu- 
lation method, the model would be placed in parallel operation with 
the system. The adaptive loop output would be applied to the model 
parameters to cause the model to adapt to the system. Thus, knowledge 
of the unadapted model governing equations combined witb the magnitude 


of the adaptive loop output would serve to define the system. 


This method is shown in block diagram form in Fíg. 13. 


Fig. 13 
MODEL MODULATION WITH MODEL ADAPTATION 
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Thís method of system identification does not relieve the engineer 
of the necessity of system analysis. The form of the transfer function, 
the order of the operating equation, and tbe general magnitude of 
the poles and zeros of the system must be known in order to synthesize 


OR 





the model. Ideally, one or two critical parameters whose magnitude 
may not be susceptible to analysis may be estimated, then measured 
or identified by tbe model modulation. These limitations should be 
expected from any adaptive system. The knowledge required to build 
an adaptive loop which would assume no knowledge of the system, and 
then synthesize perfect identification, with sufficient economy to 


be practical, is beyond the state of the art. 


The standard second order servo or feedback system equation is 


formulated in transfer function terms as: 


1 
Со) АС “Wy 


Z —-——Y 3. 
1% (9) 52 + Alu, + 27 


For this system the standard open loop transfer function is: 


Cua - G = ER 
መመመ < о) _ _ QA — 3.2 
€ | 5(5+ ደያጭነ 


The response of the second order system is completely described by 

the two variables, zeta and ሠ A higher order system, or any 

system with closed loop transfer function zeros, naturally, has 
additional significant variables. If, however, it is considered that 
there are two variables of primary interest, the higher order open loop 


transfer function may be written as: 


ኦላ 
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where ш and zeta are these variables. 
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In the development of model modulation system identification method 
two major assumptions were made. The first of these is, 
1. 
That in the operating system to be identified there will be but 
two unknown parameters. Thus the system must be second order, 
a completely dominant second order approximation of a higher 
order system, or a higher order system which has the cther 
parameters fully identified Jf the latter is the case, the 
format of the preceding equation will be used to describe the 
system. 
This first assumption limited the parameter adaptation problem to the 


two variables zeta and E 


Examination of any physical system reveals that zeta and us may 
or may not be functionally independent The simple second order system 
detailed in Fig. 5 serves as an adequate example. The closed loop 


transfer function is: 


C 4” 
5,2 % 
a ИН” 3.3 
TUS. c^ ተ T , HK — М Ка Я 
т Eu 


In this system variation of the friction(t) will yield an independent 
zeta and DM. Variation of all other factors will result in zeta and 
ሠ being different functions of the same variable, and thus not in- 
dependent. Thus zeta and » may be independent or both different 


functions of the same variable. 


30. 











The adaptive method that has been discussed to this point has 
dealt with but one variable. Since the interdependence of zeta and w 
n 


is at the most one of their being two different functions of one factor, 


it is reasonable to expect that two adaptive loops will be required. 


The second major assumption required depends heuristically on this 
discussion of the functional relationship of zeta and о This 
assumption is: 

Ji. 

That model may be adapted to the system through two separate 

adaptive loops. The first of the loops will operate upon the 

System and model output to develop the first corrective factor. 

The second adaptive loop will operate upon the derivative with 

respect to time of the system and model output to develop the 

Second corrective factor. 

The remainder of this report is devoted basically to showing that these 


two assumptions are justified. 


Proceeding on these assumptions the investigation utilized models 


of the following general form: 
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A schematic of the operation of these loops as tentatively applied to 


system identification is illustrated in Fig. 14. 
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Fig. 14 
SCHEMATIC OF TWO PARAMETER CORRECTION 
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With this model the characteristic equation of the operating 
system could then be easily detailed. For the model with the form 


of equation 3.5 and unity feedback, the closed loop transfer function 


becomes: 


Сы . Кбе)? (5а) 
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This function may then be used to accurately detail the system differ- 


ential equations. 


Only limited single parameter, single adaptive loop correction 
investigation was made. Primary investigation concerned the two par- 


ameters of the previous equation as in Fig. 14. 


Fig. 14 also illustrates the visualized mechanization of the model 


modulation identification system. Potentially the model and the adaptive 


E 





loops would be constructed through use of standard operational amplifier 
networks. Thus any moderate sized analog computer operating in con- 
junction with the scaled input and output of a physical system would 


serve as an identifier. 


idc 


Chapter IV 


ANALOG COMPUTER INVESTIGATION 


Model Modulation identification of an Equation of Known Order 





The primary investigation was centered around use of the analog 
computer. The investigation concerned use of two undetermined par- 
ameters in the system. The simulated control system transfer function 
was of the form shown in formula 3.3 The mode] transfer function ас 
of the form shown in formula 3.4. In both cases the actual computer 
arrangement was that of the transfer function format. The system 
values of VR and zeta were set up through switching circuits so that 
three different values of each could be used directly. The adaptive 
loops used on the computer were those shown in Fig. 6 In general, 
both filters in the adaptive loops were set for the break at 
omega = .1 radian. Fig. 15 presents the full analog computer diagram 
for a simulated third order system with a model tracking two unknown 
parameters. This diagram conforms to the format of formula 3.3 and 3.5. 
The initial investigation with this computer arrangement concerned a 


control system with simulated open loop transfer function of: 
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and a model transfer function: 
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The input was primarily deterministic. It consisted of a square 
wave with a frequency which was less than the bandwidth. Tests were 
also carried out however, using a random input. This random input 
was the output of a Gaussian "white" noise generator with an 18 db 
cut off at approximately the system bandwidth. At the completion of 
each run, the recorded m. and m. Were used in conjunction with the 
known model characteristics to derive the closed loop characteristic 
equation of the simulated system. This allowed comparison of actual 
System zeta and ш, with those derived from that identified as the 
system. A non-dimensional corrective parameter was examined and 
used as a partial criterion of identification. This corrective 


parameter a is best understood through examination of Fig. 16. 


Fig. 16 
CORRECTIVE PARAMETER 
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la CRIS figure, (1) is the dominant root location of the simulated 
control system, (2) is the root location of the uncorrected model, 
and (3) is the position of the corrected model root. The arithmetic 
mean à and the standard deviation, O of a, were computed for each 


Series of tests. 
РЕ 36- 





The data for the first series of tests isshown in Table I. In 
this table the variance of the tracked parameters is shown. The 
roots, zeta and ሠ have been shown for only the two extreme cases 
in each test, that of the maximum and minimum zeta. The correction 
parameter a and the sigma have been calculated from the assembled data 
of the eight situations for each test where tbe parameters do not 
match. Fig. 17 shows a sample of the analog computer output. This 
data is that of the first and last cases tabulated in Table 1 
Figure 17D will be used to illustrate the method of computing the 
characteristic equation roots of the operating control system. For 


this example, the model closed loop transfer function was 


10.7 (5+15) (10-ጠ ) 


6:53. н оиа ЯВ 
R m s (s+15) (s+104m . ) + 10.7 (s+14) (10-m_) ^ 


From figure 17D the values of the correction factors may be seen to 


be ne 1.5 and m. = 1.080. The characteristic equation is then 


5? + 26.08 s + 289.1 8 +1724. — 


The roots of this equation were obtained through digital computer 


service. These roots are: 


-5.894 + 3 9.268 
-14.29 


Fig. 18 shows a plot of the characteristic equation dominant root 
positions for the actual system and the identified system detailed 


PHUSEIg. JT. 


— 








Equation Order Approxima ted 


During the process of analyzing an operating system, it is 
inevitable that approximations would be made. These approximations 
will generally mask or ignore the very small time constants, which 
іп the s plane are the poles of large magnitude. In trying to use 
the model modulation method to identify an operating system it is 
then inevitable that the model and the system would be of different 
order. For a step input, the effect of these poles would be greatest 
during the initial portion of the transient response, and the effect 
would become inconsequential in the later portions of the response. 
It was hypothesized that if the uncertain parameters of the operating 
System concerned dominant roots, satisfactory identification would 
take place. The restriction to this hypothesis was two fold. It 
should be expected that if the input was of a random noise nature 
with no stead state, no identification would be possible unless, as 
in the first test of Table I, the equation order and the additional 
parameters are known. Additionally, with deterministic inputs, 
identification could take place only when these unknown poles were 
of sufficient magnitude to be unimportant with respect to the dominant 


response. 


To examine this restriction the tests detailed in Table II were 
hold. Since the systems to be tested were set up in transfer function 
format on the analog computer, each pole was a new block set into the 
basic arrangement of Fig. 14. The information in Table II is of the 
Same nature as in Table I, but as all nine system configurations differed 
from those of the model, all nine situations were used to calculate 
the a and C. 
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The results of these tests were as expected The random input 
revealed, as expected, lack of identification. Under deterministic 
inputs, when the postulated restrictions were held, satisfactory 
identification took place. One interesting point revealed in this 
investigation was that an adequate test for satisfactory identification 
would be made if the frequency of the input could be varied. Naturally 
this violates the basic restriction that led to the development of 
this method, i.e., no violation of the natural operation of the system, 
but is a point of interest. For this process the frequency was varied 
between a lower frequency with a period equal to the second order 
time constant and an upper frequency equal to the bandwidth. When 
the adaptive loop outputs remained relatively constant through this 
frequency range it was found that these outputs served to adequately 
identify the dominant roots when the characteristic roots were found. 
If, however, the identification varied through this range, the effect 


of the ignored poles precluded satisfactory identifícation. 


Figure 19 presents íllustrative data from the first tablulated 
case of Table II for this series of tests. Fig. 20 shows the dominant 
roots of the system for each of the tests illustrated in Fig. 19, 


with the dominant roots as identified from the data. 


Examination of the data reveals that for both of these types ot 
equation situations good dominant root identification has been achieved. 
It was found that these dominant roots could be identified to within 
an average error of 2.2% of the undamped natural frequency and to 


within an average error of 4.76% of the complex phase angle. When 
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this method is used to identify the characteristics of an actual linear 
control system, it is highly probable that greater accuracy may be 
attained since the model parameters may be varied to obtain error 


bracketing. This process is illustrated in Chapter VI. 
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Chapter V 


DIGITAL COMPUTER SIMULATION 


Direct utilization of a digital computer for identification of 
an operating control system was not considered It is possible, 
however, that operating systems might be monitored or controlled by 
digital techniques. In this case the output would be presented in 
the form of discrete numerical data. Ideally this data could be used 
as input to a digital-to-analog converter, and the analog techniques 
previously described could be used for operating system identification. 
In the absence of this equipment, digital simulation of the control 
system, model, and adaptive loops was investigated to illustrate the 
techniques and the potential applicability of model modulation techniques 
through numerical techniques. This investigation was conducted through 
the use of the CDC 1604 Digital Computer. The program code was written 
in Fortran. A Standard Runge-Kutta numerical integration process was 
used to solve the differential equation which governed the operation 
of the full model, system and adaptive loop. The differential equations 


were formulated in a transfer function format. 


Three basic configurations of model adaptation were tested through 
this digital computer simulation. For these configurations the model 
was used to identify the operating control system characteristic 
equation roots under the following conditions: 

1. Second order control system with undetermined pole. 


2. Second order model, third order control system with uncertain 
gain and uncertain dominant pole. 


3. Third order system and model with two uncertain poles. 





The actual investigation was limited to achieving a successful solution 
for one parameter set for each of these configurations. The excessive 
time used by integration techniques on the digital computer was of 
course a limiting factor. The first of these configurations will be 
by-passed in this discussion since it can be considered to be a sub- 


system of the second configuration. 


The model for the second configuration had the following open 
loop transfer function; 


10 (10 - m.) 


— s (s + 10 + m.) 3.1 


and a closed loop transfer function of 


10 (I0 -m ) 
A MEE ም 5.2 
R "m s (s * 10 + m.) +10 (10 = m.) | 
The control system open loop transfer functions were 
12 ООО 
Gis = вее = Dm 
s (s+100)(s+8) 
and 
Gu) = _ Bs — 5 4 
5 OOS 
Figure 21 presents the transfer function format which was used to 
formulate the differential equation for the FORTRAN program. This 
format allowed formulation of the operation by a set of quasi-state 
variables of the form: 
x = f (хү Ха» X3> 000 X ? R) O 
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The solution of the equations was programmed to be presented in both 
numerical and graphical format. Fig. 22 presents a sample of the 

input program Fig. 23 presents a sample of the numerical data output, 
and Fig. 24 presents the graphical data output for the control system 


open loop transfer function given by equation 5 3. 


The graph of Fig. 24 presents the corrective parameters as 


"corrected pole" = (10 + m, ) and "corrected gain" = (10 - т.) 
Table III shows the identification of the characteristic equation 


roots for these two situations. 


Table IIIA 
DIGITAL COMPUTER CHARACTERISTIC EQUATION ROOTS 


Control System Modél as- Corrected 
Edüation 5.3 = omo queas X D quy 

- 2101,27 5101919 
Equation 5.4 "d DS mI DD "8-00 E 

- 100.9 - 100.93 
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The third configuration examined was that for the following model 


open loop transfer function: 


G(s) » 10,000 5 
s (s + 100 + m.) (s + 10 + me) | 


with control system open loop transfer functions of: 


10,000 
И а р ааа 8 s; E: 
ፅ እ ም ጌን ን ቢቢ ዓዓ ыт) 2 
and 
ды 10,000 з 


Sis +t OJ (STE 725) 


The equations were formulated in the manner of Fig. 20, except 

that the adaptive loops each utilized one filter and one integrator 
instead of the two filters previously detailed. Satisfactory results 
were obtained. Fig. 25 presents sample results from the investigation 
of equation 5.7, and Table IV presents the root identification 


achieved. 


Table IIIB 
DIGITAL COMPUTATION CHARACTERISTIC EQUATION ROOTS 
TWO UNCERTAIN POLES 


Control System Model as Corrected 
Equation 5.7 smog 25 sv 5.72 ISI ENE 

- 150,49 -150.49 
Equation 5.8 = ሪር) ት. 11.1) 290085 + j 11.05 

- 76.87 - 74. 56 
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Chapter VI 


ACTUAL SERVOMECHANISM IDENTIFICATION 


Introduction 
The previous discussion and the data in Table I and Table II 
indicate that the identified location of the roots of an operating 
system is the function of many variables. The type of model and 
the proximity of the model equation roots to the rcots of the equation 
of the operating system are obviously factors of importance in this 
identification. With this point in mind, it is suggested that the 
following steps be followed in the process of identifying an actual 
system. 
l. Analyze system operating factors to enable determination of 
open and closed loop transfer functions. 
2. Determine which values in the open loop transfer function 
are questionable. 
3. Utilize operational amplifiers or an analog computer to 
prepare the dynamic model and the adaptive loops. 
4 Scale the model input and outputs to the estimated 
corresponding values of the system. 
5. Set indicated adaptive loop polarities to adapt the model 
to the undetermined or questionable values in the system. 
6. Adjust gain of adaptive loop to obtain satisfactory operation 
7. Calculate the transfer function of loop corrected model, 
scale a new model to these functions and repeat operation 
to generate better conformance or to confirm values pre- 


viously attained. 


de 





A position servomechanism was used to illustrate the applicability 
of the model modulation techniques to an actual operating system. 
The servo used was one developed for use in instructional laboratories 
by Professor J. R. Ward of the U. S. Naval Postgraduate School, and 
is partially described in Ref. 6. The full description and details 


of the operation of this system are on file at this school. 


Fig. 26 presents a block diagram of the system. This repre- 
sentation shows that the servomechanism may be considered as a gear- 
head system. Thus the use of this identification method conforms to 


the stated requirement of parallel model operation. 


Fig. 26 
POSITION SERVOMECHANISM 
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Analysis of this mechanism resulted in the transfer function diagram 


Of Fie. 2/. 


13% 





pis. 27 
POSITION SERVOMECHANISM BLOCK DIAGRAM 








The approximate parameter values in the system were determined. The 
magnitude of inductance (L) was such that the approximation L ^O 


was made. The open loop transfer function was then: 


EO 
Es RJ ረ 
— = (5 З ню о ran = оо ae 6 1 
€ RE + Kay 5 (ፄ፥ ኮ) 
(++ Rhe] й 


The value of friction(f) was difficult to evaluate because of the 
expected non-linearity of friction in a d.c. motor. A nominal value 
was chosen to represent the viscous friction, while coulomb friction 
effects were ignored. It was further assumed that the open loop gain 
was questionable. As a first estimate it was assumed that the open 


loop transfer function was: 


e = [ረ AOO 
im 211171. ^ s(s+8) er 


г 


Fig. 28 presents the control system-model-adaptive loop block 


diagram. 


ne 
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The actual model and adaptive loop analog computer schematic was that 
detailed in Fig. 15, while the scaling factors were achieved through 
cascaded operational amplifiers to assure correct polarity. The 
servomechanism and the model were driven by a square wave of one cycle 
per second. Table IV presents the results of four tests made first 
to bracket, then to closely approximate the open loop transfer 
function. Fig. 29 presents the results of the final run listed in 
Table IV. One the basis of these results the transfer function of 
this servomechanism was estimated to be: 
625 

ёс 5 + 04656 + 625 
These results compare favorably with information derived from other 


sources. Calculated results derived for careful component analysis 
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Table IV 


POSITION SERVOMECHANISM CHARACTERISTIC EQUATION IDENTIFICATION 


OPEN LOOP TRANSFER FUNCTION 


Model 


10 (40-m_) 
S (s+8+m.) 


10 (70-m_) 


s (s+10+m,) 


10 (55-m ) 


10 (60-m ) 


S (stlOtm.) 


S (st10.15) 


SYSTEM DYNAMIC EQUATION 


Corrected Model Identified Roots 
574 -4,83 +3 23.45 
s (s+9.66) 


С = ‚201 ш = 24 
n 


652 -5.075 +j 25 
ር = .204 9 = 25.55 


622 E5517 !፡!:. 22 


s (st10.24) ር = .206 9 - 24.95 
n 


625 -5.075 +3 24.4 


s (s+10.15) C= .208 ወ = 25 
| п 
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indicated a closed loop transfer function of 


— 


e, (s) 590 
е; 5) 555 (12.5 %1.5)6 + 590 





while various open loop frequency response tests indicated the closed 


loop response to be: 


Ce (s) ES 5695 f 70 
e, (s) s? 4 (10 t 1)s + 565% 70 


6.5 


In view of the non-linearities involved in a d.c. servo motor it was 
felt that these model identified roots represented a successful 


illustration of the validity of the method. 


Discussion 

These servomechanism tests confirmed that the model modulation 
method could be successfully applied to obtain the governing operating 
equations of a servomechanism. Beyond this the tests emphasized the 
necessity of continuous excitation of the control system for successful 
model adaptation. For a linear control system the model modulation 
adaptation will operate only while the servo is undergoing excitated 
or active response to an input signal. It would be expected that the 
correction or adaptive signal would be generated in response to a 
step for a period equal to the settling time of the system. Beyond 
the settling time, with no excitation response operation to generate 
error, the time constants of the adaptive filter network will allow 
dissipation of the adaptive signal magnitude. The analog computer 


simulations detailed in this chapter illustrate this operation of the 
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adaptive loop For the analog tests the general input frequency of 
the square wave was 0.5 cps. The average settling time for the 
second order model was T u 1.0 sec, and the bandwidth was 

f ™ 1.9 cps. When tests were run with systems with small zeta, lower 
input frequencies would achieve successful adaptation. In all cases, 
successful operation was attained with frequencies which were as high 


as the bandwidth frequency. 


In the case of the operating servo, the non-linearities of the 
friction caused a faster damping than the linear interpolation would 
indicate. For this system, the adaptation proceeded successfully 
when the driving square wave frequency allowed limited oscillation 
per step. This in effect was a partial masking of the non-linearity 


through testing techniques. 


These situations emphasize the limitation of control system 
excitation in the use of the model modulation method. It is obvious 
that this limitation exists for both system adaptation and for system 
identification through model adaptation. A potential solution of 
this problem might be found through investigation into utilization 
of signal identification circuits. This signal identification might 
be applied to the output of the control system to control the operation 
of the adaptive circuit. This identifier could allow the adaptive 
circuit to operate while the control system would be underdoing 
oscillatory or excited response. When this control system response 
became quiescent, the identifier could place the adaptive loop into 


a hold state with a constant level of output maintained to maintain 


“ርያ 





the adaptation until such time as the control system would be again 
activated. Additionally, it is expected that the adaptive loop 

elements would be synthesized through the use of operational ampli- 
fiers. Since these are magnitude sensitive, another role of the signal 
identifier would be to hold adaptive loop operation when signal  mag- 


nitude became excessive. 


The techniques of sampled-data control systems offer attractive 
possibilities for circuits to fulfill these requirements. Hybrid 
digital-analog circuits might also be investigated. In any case, 


these techniques warrant additional investigation. 
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Chapter VII 


CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 

It is concluded from the results of this investigation that the 
identification of the characteristic equation roots of a control 
system by model adaptation through model modulation is a valid 
technique. The technique is subject to operating restrictions as 
follows: 

1. Тһе control system must be linear or be subjected to operating 
non-linearities of minor magnitude. 

2. The open loop transfer function must have but one or two 
parameters whose magnitude is uncertain. 

3. The additional parameters of the open loop transfer function 
must be known or be of sufficient magnitude to be incon- 
sequential. 

4. The control system must be undergoing active operation in 


response to an input signal. 


It is further concluded that: 

1. The roots of the characteristic equation of a control system 
may be identified with a high degree of accuracy. During 
these tests, these dominant roots were identified to within 
an average error of 2.2% of the undamped natural frequency 
and to within an average error of 4.76% of the complex phase 


angle. 


а 





2. When the order of the characteristic equation is unknown, 
second order dominant mode approximations will yield 


consistent results which are valid in the approximation 


Recommendations 
The results of this investigation indicate that there are 

certain areas of interest which warrant further investigation. The 
first of these is investigation into the stability of the dynamic 
adaptive loop. The stability limits are subject to various major 
bounds in loop gain, input signal frequency, and input signal magnitude. 
Theoretical analysis of these bounds will be extremely difficult, 
thus analog computer investigation might offer a valid method of 


investigation. 


The second of these would be the investigation of methods to 
remove the necessity of a continual excited response from the control 
system. Certain suggestions have been made in the text which might 


be applied to achieve this end. 


The third suggested area of investigation is in the identifi- 


cation of control system operation subject to known non-linearities. 


— 
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